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Abstract 

The boundary problem of linear classical optics about the interaction of electromag- 
netic radiation with a thin dielectric film has been solved under explicit consideration of its 
discrete structure. The main attention has been paid to the investigation of the near-zone 
optical response of dielectrics. The laws of reflection and refraction for discrete structures 
CN ■ in the case of a regular atomic distribution are studied and the structure of evanescent 

harmonics induced by an external plane wave near the surface is investigated in details. 
It is shown by means of analytical and numerical calculations that due to the existence of 
the evanescent harmonics the laws of reflection and refraction at the distances from the 
surface less than two interatomic distances are principally different from the Presnel laws. 
From the practical point of view the results of this work might be useful for the near-field 
^h * 1 optical microscopy of ultrahigh resolution. 



PACS numbers: 78.20.-e, 81.40.Tv, 42.25.-p, 87.64.Xx 



; 1 Introduction 

Usually in the investigations of optical phenomena in dielectrics a medium is considered as 
a continuous system of radiators. Besides the macroscopic Maxwell equations are employed. 
This way of description works well if one considers the phenomena which take place at large 
enough distances from medium interfaces. However, in recent years the interest to the optical 
phenomena near the surfaces of different media at the distances much less than the radiation 
wavelength has dramatically increased. Beside a purely fundamental interest this is connected 
to the vigorous development of the scanning near-field optical microscopy 0-0. There were 
developed various schemes of optical near-field microscopes and the work in this direction is still 
currently underway. However, all the schemes have a common feature, namely, the measurement 
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of the optical response of the medium on the external radiation is carried out in the near zone 
at the distances of the order of surface inhomogenity j|, 0] . 

The resolution obtained with the aid of near-field optical microscopes permanently increases 
and at the moment it is of the order of 10 nm. Besides the probe scans the surface also at 
the distance of the order of 10 nm 0. At such a small distance it becomes necessary to take 
into account the discrete structure of the medium. The attempts of consideration of discrete 
medium structure in the theoretical description of the near-field optical measurements have been 
undertaken in recent years ||. However, the main attention has been paid to the situations 
when discrete dipoles model rather large surface inhomogenities and the spatial distribution of 
these dipoles is asymmetric. The progress in the near-field optical microscopy allows one to hope 
that there will be possible to observe the atomic structure of surfaces. Then it will be necessary 
to interpret the optical response from single atoms more or less regularly distributed on the 
surfaces. Similar problems occur when one wants to describe the near-field in the vicinity of 
the surfaces of photon || and plasmon || crystals. Besides it becomes more natural to employ 
microscopic equations. 

Intensive studies on optics of discrete dielectric media have been started not very long time 
ago. The laws of refraction and reflection for discrete structures were shown to be different 
from Fresnel laws [JT0|] and this leads to the anomalies in the behavior of the electromagnetic 
field near the surface |TT|-|jr3[. Theorem of extinction for discrete structures has been derived 
and the corrections to the Snell law were worked out HI4" 
spontaneous emission process [15|, |l( 



The characteristic features of the 



structures were investigated 
considered in Ref. |19| 



, superradiation [|l7j and optical bistability [|18[] in discrete 
The second harmonic generation in discrete structures has been 
where a good agreement with experimental data has been pointed out. 



It was also shown that the consideration of the discrete structure of the medium leads to the 
surface local-field effect [^(J and surface induced optical anisotropy pl|-p3 . 

Properties of the waves reflected from discrete structures in the wave and near zones have 
been studied as well. It has been shown that in the wave zone the reflected waves can be 

the treatment in Refs. 



described by plane waves [JlOj, [LI], |25|, |26], ^7 



25|, |26 



does 



However, 

not take into account evanescent harmonics which give a considerable contribution into optical 
response of discrete dielectric structures in the near zone. These evanescent harmonics are not 
related to the total internal reflection and they exist even in the case of a single atom P5 |. Due 
to the neglect of the evanescent harmonics the method developed in Refs. p5| , 26 , ^% can be 



employed for the description of the optical response of a thick enough medium in the wave zone 
only. 

A more general approach which implicitly takes into account the evanescent harmonics has 
been developed in Refs. |L0], [11]. This allowed us to show that the fields of the reflected waves 
in the near zone are periodic functions of longitudinal coordinates with the period equal to 
the lattice constant. This result confirms the assumption put forward in Ref. that with 
the aid of near-field optical microscopes one can observe the atomic structure of the medium 
surface. In the present paper we shall continue the investigations started in Refs. 
Employing the methods of the calculation of the lattice sums developed in Refs. 



10, 11 



12, 24 



we 

are going to investigate in more details the behavior of evanescent harmonics near the surface 
of a dielectric film with an ideal crystal structure. We would like to note that the role of the 
evanescent harmonics in the near-zone optical response of dielectrics has been point out by 
Sivukhin |29fl. On the basis of rather simple and intuitive physical arguments he came to the 
qualitative conclusion that the Fresnel laws of reflection and refraction have to be violated at 
the distances from the surface less than one interatomic distance. In the present paper we are 
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going to undertake a quantitative analysis of this fundamental issue. Besides it will be shown 
that the existence of the evanescent harmonics leads to principal deviations from Fresnel laws 
of reflection and refraction at the distances from the surface of the order of few interatomic 
distances. 

The paper is organized as follows. In section 2 we shall discuss the problem we are going 
to solve and the main equations which constitute the basis of the microscopic approach to the 
description of optical phenomena in the linear classical optics. Section 3 is devoted to the 
analyses of the probe influence on the electromagnetic field distribution in the film and the 
back influence of the perturbation caused by the probe on the field magnitude on it. It is 
shown that in some cases the probe influence can be neglected. In section 4 we shall calculate 
the distribution of the local field in a dielectric film under the incidence of a plane wave. In 
section 5 we shall analyze in details the characteristic features in the behavior of reflected and 
transmitted waves in the wave and near zones. 



2 Statement of the problem and basic equations 

Let a monochromatic light wave with frequency lu and electric field strength vector E/(r,i) = 
Aj(r) exp (— iut) be incident on a dielectric film with thickness h. Above the film there is a 
probe which measures the electromagnetic field at a certain point of space. The probe will be 
treated as a single point-like dipole with a linear polarizability a p . Such a scheme can be realized 
in a real experiment using either a single atom in a magneto-optical trap [^] or an admixed 
atom implanted in the needle of a near-field microscope |7|, |31]| . Our aim is to investigate 
the behavior of the electromagnetic field in the film and to calculate the characteristics of the 
electromagnetic field on the probe that is moved (scanned) along the film surface at some fixed 
distance d. The field detection, which represents itself a separate problem, will not be touched 
in the present paper. 

We shall treat the dielectric film as a system of N monolayers located at a distance a from one 
another. The atoms inside the monolayers possess a discrete spatial distribution. In a stationary 
case the equation for the local electric field strength on the probe Ep(r, t) = E p (r)exp (— iuit) 



32], |33| can be written down in the form 

2 N 

E p (r) = Aj(r) + t-kla p E p (r) + a ^^V r x V r x E^rf.)^), (1) 

where r = (x, y, z) is a radius- vector of the probe position; a is the microscopic polarizabil- 
ity of the atoms inside the film, which is assumed to be independent of the field. G(R) = 
exp (ikoR) /R is the Green function of Helmholtz equation, k = uj/c is the vacuum wave 
number, H aj = r — v ap R a . = H aj , r a . = (rjj., Zj). We have also introduced the notation 



E,.(r« .) = E(r» ., Zj ), Zj = -(j - l)a, j = 1, N. (2) 



The reference frame is chosen in such a manner that z axis is perpendicular to the film surface. 
The second term on the r.h.s. of eq.([l|) describes the Lorentz radiation damping [32, Q. 



In spite of the fact that the radiation damping makes usually a small contribution into the 
resultant field the consideration of this term is of principal significance. If one does not take it 
into account the energy conservation is violated (see []35| and references therein). 
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Let r„, — (x an y ai , z\) be a radius- vector of some atom in the Z-th monolayer. Then for the 
local field at the position of this atom we have 

E ; ( r y = E n (rjlj + ihlaE t (rjj + a p V r<l; x V ra( x E p {v)G(R ai ) 

N 

+ «EEV raj x V ra; x Ej(rl.)G(R aia .), l=TJ7, (3) 
3=1 

Note that the term with r Qj = r a[ has to be 
excluded from the summation. 

We can look for the solution in the l-th monolayer of the film in the form 

^( r l)=E?(rl)+Ef(ri), (4) 
and in addition we impose a requirement that the field E° (rj| ; ) has to satisfy the equation 

2 N 

E° (rjl ) = E n (rj,) + «-fc 3 aE? (rj,) + « £ £ Vr a; x V ra; x EjCrJ^GCi^o,). (5) 

6 3=1 a j 

Then from eqs.(|l|), @ it follows that the fields Ef and E p will be the solutions of the system 
of equations 

2 N 
E p (r)=A 7 (r)+^-/c 3 a p E p (r) + «^^V r x V r x E»(r» 

3=1 a o 
N 

+ «EEV r xV r xEj(rl)G(iJ % ), (6) 

3=1 a i 

a P V rai x V rai x E p (r)G( J R a J 

v 

«EEVr, x V ro; x E^)G(R aiaj ). (7) 

3=1 a J 

The physical meaning of the fields E[* and Ef is the following. E° is the local field inside the 
film without consideration of the contribution made by the probe. This contribution is taken 
into account by means of the field Ef. Thus, for the solution of the problem stated above it is 
necessary first to solve the unperturbed problem @. Then using the obtained solution one can 
determine the perturbation of the field inside the film caused by the probe and calculate the 
back influence of this perturbation on the field at the position of the probe (f|), (0). However, 
before dealing with this general problem we shall estimate the perturbation in the film caused 
by the probe and find out under which conditions it is considerable. 

3 Estimate of the perturbation caused by the probe 

In this section we shall estimate Ef which enters eqs.(|]), (0). For that purpose one can neglect 
radiation damping terms in eqs.(||), (0). Our starting point is the following. The perturbation 



where En (r|[J = A 7 (r a; ), R aiaj 
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of the field in the medium Ef , which is caused by the probe, decreases rapidly with the increase 
of the distance between the surface and the probe. If in addition we take into account that the 
influence of atoms on one another also rapidly decreases with the distance, one can come to the 
conclusion that either the probe or the nearby atoms will be influenced mainly by the atoms 
located in the nearest vicinity around them. Thus, it is necessary to find the field E p (r a ) in 
the direct vicinity of the probe. Let's estimate first the size of this region. 

One can assume that the behavior of the field Ef has the same features as that induced by 
the probe 



Ef(r«) 



oc 



3(E p R ai/ 



R 



a, 




R 5 

EpR<j; )R-a; 



R 3 




zk ( 3 - {EpKa > 



)R 



•Oj 



i? 4 



(8) 



— — exp(ik R, 



ai, 



Taking into account that the field on the probe is defined by the atoms in the near zone, one can 
keep in eq.([8]) only the terms proportional to -R^ 3 - We substitute the corresponding expression 
into (|7|) and replace the summation by integration. As a result of this the expression for the 
third term in the r.h.s. of ([?]) takes the form 



N 2 

« E E V ra; x V ra; x mv\)G{R aiaj ) -+p[ V r x V r x EV)dr', 

3=1 a J Ect 



(9) 



with p being the macroscopic concentration of atoms. The integration in ([|) is carried out over 
the film volume limited by the external surface S. In order to prevent the interaction of the 
atom with itself a small region limited by the spherical surface £ CT centered at the observation 
point has to be excluded from the integration. The size a of this small region is of the order of 



lattice constant. In Ref. |10[ it was shown that in order to achieve the complete correspondence 



f a. 



passing from discrete to continuous atomic distribution one has to put a 

Let's determine the region of space, atoms from which make the main contribution to the 
field on the atom nearest to to the probe. Let 5 be a tolerance of the calculations of the dipole 
field. Then 



1 - 



V r x V r x E p (r')dr' 




x V r x W(r')dr' = 5, 



(10) 



Integration in the denominator and in the numerator in eq.(|l^) is carried out over the film 
volume and over the spherical segment with the radius L centered at the position of the atom 
nearest to the probe, respectively. Taking into account (§) we obtain from eq.flTOD the following 
estimate for the radius of the sphere L within which it is necessary to calculate the fields on 
the atoms: 

L « (Tt 1/3 . (11) 

Let 5 = 0.1%, then L w 7a. Note that in the calculations one can use even less values of 
L, because in eq.(^) there is also the field induced by the probe and compared to this field 
the influence of the remote atoms is less. The conclusion that only the neighboring atoms 
contribute to the field at the site of probe in near-field regime is in agreement with the results 
of Ref. 



36 



where the influence on the probe caused by microscopic dielectric spheres located 
at the surface has been calculated. It was shown that the lattice of 5 x 5 spheres makes the 
main contribution. The increase of this area does not influence the final result. 



5 



Performing analogous calculations for eq.@ we obtain that the main contribution to the 
field on the probe is made by the atoms located inside the sphere with the radius L centered 
at the probe position. Thus, the system of equations (||), (0) is reduced to the system of few 
linear algebraic equations, which can be solved numerically after the field E° is obtained from 
eq.(|5]). However, as it will be shown below in many situations it is even not necessary to solve 
this system of equations due to the negligible contribution of E p . 

For the sake of distinctness we assume the probe to be located above some surface atom 
of the film. The field induced by the probe at the position of this atom is determined by the 
expression 

E e// « -a p E p /d 3 , (12) 



where d is the distance between the probe and the atom. In formula fll2|) we have neglected 
the retardation and the field on the probe is chosen to be polarized parallel to the film surface 
assuming that the variations of the direction of the polarization vector does not influence much 
the final results. It has been shown above that in the calculations of the field on the probe it 
is enough to take into account the atoms located inside the sphere with the radius L. The rest 
of the film does not influence much the probe. So we shall estimate the influence on the probe, 
caused by the atoms located within the spherical segment centered at the position of the atom 
below the probe. Obviously, the microscopic field amplitude on the atom under consideration 
is larger than that on the other atoms of the film. The upper estimate of the field amplitude 
at the position of this atom E max is given by 

|E maa; | < |E e// | [1 - l.nCHL/a)]- 1 , (13) 

where the dimensionless parameter C = ap is a volumic polarizability [p7| . We shall assume 
that the field on any atom inside the spherical segment is constant and equal to the field on the 
central atom E max . It is obvious that in this approximation the influence on the probe caused 
by the spherical segment is larger than it is. Our calculations are approximate, so we replace 
the summation by the integration. Calculating the field at the position of the probe we get an 
upper estimate of the last term in the r.h.s. of eq.(H) 

|E^| < 1.71C\E max \ In (1 + L/d) . (14) 

Taking into account (0), (|13|) , (|14D, we get a condition under which one can neglect this term: 

o9_ 1.71Cln(l + L/d) 

a~ p ^ [1- 1.71Cln(L/a)]' ( } 

Note that the condition (pj|) does not depend on the way of excitation of the surface atoms. It 
remains the same if the atoms are excited by a wave incident from the side of the probe, by an 
evanescent wave, or by the wave radiated by the probe itself. 

In eq.(|HJ) ln(L/<i+ 1) and \n(L/a) are of the order of unity. Therefore, one can rewrite (|T5|) 
in a simpler form 

d 3 1.71C , , 

» ; ; 7 , r - 16 
a p 1 — 1.716 

Usually for dielectric media C ~ 0.1 and d 3 /a p ~ 10 at the distances from the surface equal 
to one lattice constant. Therefore, inequality ( |ISD is satisfied at the distances equal or greater 
than one lattice constant. 
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At the distances between the probe and the surface such that the influence of discrete atomic 
distribution is negligible but the retardation is still important the calculation of the field on the 
probe reduces to the electrostatic problem of the calculation of the dipole field near a dielectric 
surface. For the field ~Ef" at the probe position we get 

E ^TTI^ (E >- 3(E > n)n) ' (17) 

where n = (0, 0, 1) is a unit vector perpendicular to the monolayer. We write here because 
in the derivation of eq.fll7|) we have used the Lorentz-Lorenz formula for the refractive index n: 

1 + (8n/3)C 

U ~ 1 - (4vr/3)C ' (18) 



which does not describe properly optical properties in the surface region [37]. The condition 
which allows to neglect the last term in eq.(0) is given now by 

d 3 n C 

a p 4 1 + |7rC v ' 

The condition (|19"D is in a good agreement with (fL6|), albeit the former imposes less strict 
limitations on d. 

Thus, at the distances between the probe and the surface of the order of one lattice constant 
the back influence of the film on the probe can be safely neglected. In this approximation the 
field on the probe is determined by the field of the external wave and the field E°. The problem 
of calculation of the field E° will be discussed in the next section. 

4 Distribution of the field inside the film (solution of 
the unperturbed problem) 

Let the external wave be a plane one, i.e., 

A/(r) = E / exp(ik r), k = /co(sin 9/ cos sin 9/ sin $/, — cos 9/) , (20) 

where 9/ is an incident angle and the azimuthal angle defines the orientation of the incidence 
plane. 

The atoms of the monolayers are assumed to form a regular periodic structure with elemen- 
tary translations vectors a x and a 2 . Without loss of generality one can put 

ai = an (1, 0, 0), a 2 = a[|(a,/3,0). (21) 

For the sake of simplicity the vectors &i and a 2 are assumed to be the same for all monolayers. 
The generalization on different ai and a 2 for different monolayers can be easily done. In these 
notations the radius-vector of some atom in the j-th monolayer has the form 

r a . =r° j + a mn = r° + max + na 2 , (22) 

where m and n are integer. The vector has the following form 



r 



o 



ctij-ai + a 2j a 2 + zjn , < «i(2)j < 1 • (23) 



7 



Nonvanishing ai\j and a 2 j take into account possible parallel shifts of the atomic planes relative 
to one another. 

Making use of the principle of parallel translational symmetry the solution for the local field 
in the l-th monolayer can be written down in the form 



E?( 



E? 



exp ^k r' 



(24) 



We substitute (|2T)D, (^|) into (|5[). Now it is necessary to calculate the lattice sums in eq.(|5|). 
Due to the slow convergence of lattice sums their direct calculation is a very laborious and 
in fact impracticable computational problem. In order to compute them one can use "the 
Lorentz method" according to which the atoms located near the observation point are treated 
as discretely distributed and the atoms outside a fictitious boundary are treated as continuously 
distributed. This method can be employed either for well-ordered or for random media and for 
the observation points either inside or outside of the media. This way of computation has been 
employed in various modifications in Refs. JTD|, O, EUL RH]. On the other hand, if the atomic 
distribution in the monolayers is regular, this separation of discrete and continuous regions is 
not necessary anymore. In this case there exist another way of computation of the lattice sums 
which is rapidly convergent and more convenient for the interpretation of the field behavior 
near the surface [|i~2"j 24|| . It is this latter way will be employed in the present paper. 

If the observation point is located not in the monolayer over which the summation is carried 
out (j 7^ I) the lattice sums can be converted to the form which is more appropriate for numerical 
calculations, making use of the formula for the dipole field induced by the j-th monolayer at 
some point r, which has been obtained in Ref. |24j by means of Fourier transform (Ewald's 
three-fold integral transform): 



a 



EV r xV r x E°exp (ik rjl.) G 



"J 
OO 



E A p?( r ' z i) exp(ikJlr) , z ^ 



J 3 i 



(25) 



p,q= — oo 

lai x a 2 | 



kpg X (is-pq x 



exp 



(k^-kSXr-r") 



■"pq 



= | (kjj + g| 9 , Kpq ) if z>Zj , 2g j 

Kpq = \Jkl~ (kj + gpg) 2 , gjg = pgi + qg 2 (27) 
Icq = /c ( sm ©7" cos sm @/ si n 0). 

The vectors of the reciprocal lattice gi and g 2 are related to the vectors a x and a 2 as 

a 2 x n n x ai . 
gl = 2tt^ , , g 2 = 2n- 28 

I a l x a 2| | a l x a 2| 

The choice of the vectors ai and a 2 can be arbitrary. However, it makes sense to choose them 
in such a manner that they have a minimal length. Then the length of the vectors gi and 
g 2 will be also minimal. Such a choice of the basis is convenient, because in this case the 
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decay coefficients of the evanescent waves for square lattices and for slightly anisotropic lattices 
increase with the parameter k = p 2 + q 2 . 

Formula (|25|) is a useful representation of the lattice sums for numerical calculations because 
the sums over p, q converge much faster than the sums over aj, provided that \z — Zj\ > an. 
The calculations carried out in the next section show that in the sum over p, q it is enough to 
keep only few summands with minimal g^ q . 



On the other hand, formula (25) has a very clear physical interpretation. It represents itself 
a decomposition of the field induced by the monolayer of discrete atoms on propagating and 



evanescent waves. Indeed, if 



K + 



< k the quantity n pq is real and we have propagating 



waves. Otherwise, n pq is imaginary and we have exponentially decaying waves. A more detailed 
discussion of the role of propagating and evanescent waves will be given in the next section. 

Note that the quantities A pq (r, Zj) with p, q ^ 0, i.e. the evanescent harmonics, are periodic 
functions of x and y and they posses the following property of the translational symmetry 



A p „(r + ma 1 + na 2 , z 7 - 



A pi j(r, Zj j 



(29) 



As we shall see further on, this property of the evanescent harmonics leads to the fact that the 
intensities of the reflected and transmitted waves near the film surface are periodic functions 
of x and y. 

If the observation point r is at the position of some atom inside the monolayer over which 
the summation is carried out (j = I), then for the calculation of the lattice sum it is convenient 
to employ the method based on Ewald's onefold integral transform [12]. Separating the initial 
sum over the real lattice on the sums over the real and reciprocal lattices one gets the following 



expression for the dipole field [12 



a£V r xV r x E° exp (ik rj . ) G (R a . ) 



+ z-^aE°exp(*k|jr) = a/(k )E° exp(zkj|r) , (30) 



where the term with r" = rj|. has to be excluded from the summation. Tensor / is a symmetric 



one and the components f xz , f yz , f zx , f zy vanish. It has the form |12] 



VK 



ai x a 2 



E 

p,q 



k 2 5-, 



1 exp(zka r , 
+ 2 ^ ~a? 



+r 4 

1 r. 



- k k 

A M (1 + T) 

Z P1 



7]H 



P'-l 



If 1 \fi F 2 A- t\ v a M F 

1 nm y J i 1 ^ 1 - nm ' nm nm 



pV a M 

nm nm nm 



+ C.C.} 



3 

nm 



-ik 3 



(31) 



where u, // = x, y, z, a Tl 



A pq = erfc 



»? = M«: r=(-l) 5 -(-l) 

ik, 



2E 



lAE ( k 2 pq " 



= exp (-ika nm ) erfc I a nm E - ^ 



1 ika nm + k a 
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-k 2 + 



3ik 



+ 



2E 



71 



exp 



_zp2 2 i _ 

- C/ u nm ~r 4^2 



-\- IE Q nm . 



The rate of convergence of the sums in ( |31~D is determined by the parameter E, which has a 
dimension of the inverse length. It appears as a result of the formal separation of the initial sum 
over the real lattice on the sums over the real and reciprocal lattices. This parameter has no any 



physical meaning, it must be only real and positive and its choice can be arbitrary [T2|]. However, 
it makes sense to choose the parameter E in such a manner that the rates of convergence for 
the sums over (p, q) and (m, n) are the same. One can show that this requirement is fulfilled, 
provided that E = ^Jtt/\sli x a 2 |. Substituting this value into (|3"ID one can estimate the maximal 
values of (n, m) and (p, q) which has to be taken for the calculation of the lattice sums with 
the tolerance e. Taking into account that at large x erfc(a;) ~ exp(— x 2 )/ (xy/rr), we get the 
condition 

,2 



exp^-7ra; m /|ai x a 2 |) « e. (32) 

Exact numerical calculations for a square lattice show that if the summation is restricted by 
the terms with n, m,p,q = 0, ±1 the relative tolerance of the calculations of the dipole field is 
about 10~ 6 , which is in a good agreement with the estimate ([^). Therefore, our choice of the 
parameter E appears to be very convenient for earring out numerical calculations. 

Note that in the computer calculations of the dipole field the algorithm based on Ewald's 
onefold integral transform is 5-10 times faster than that based on the Lorentz method. 

Taking into account identities 
following form: 



5]), (|30| ) the system of equations (|J) can be rewritten in the 



«/(k ) 



E 



N oo 

E n - i2nCa £ ' £ [k pq x (k pq x E° 

j = l p,q=-oo 



exp 



(k P g 



c?) 



pi 



where C = a/{a\a.\ x a 2 | 



E 



n 



E, 



exp (ik t 



T ai 



Eqj exp ( — iko cos 0/^;) , l = l,N. 



(33) 



(34) 



Prime at the sum sign in ( |3"3"D means that the summand with j = I is excluded. 

The results of numerical solution of the system of equations ([33|) for the cubic lattice and 
s-polarized incident wave are shown in Fig. |l|. In order to demonstrate to the fullest extent the 
characteristic features of the field behavior near the film surface we have used a rather high value 
of the volumic polarizability (C = 0.2). Normally in the optical range C = 0.04 0.15 which 
corresponds to the refractive index n = 1.27 -j- 2.46 ||39|| . However, for some materials C can be 
higher. For instance, na a p {^o — 564 nm) = 
which corresponds to Cc a p = 0.19 and Csi 



3.42 [0, H and n Si (X 
= 0.2, respectively. 



620 nm) = 3.94 22, 40 



In Fig. |I] the results of the calculations according to the Airy formulae are also shown for 
comparison. Besides we have used the Airy formulae for two cases: (1) the film occupies the 
region of space — a(N — 1) < z < 0, and (2) the film occupies the region of space — a(N — 1/2) < 



z < a/2. In the former case the corresponding formula has the form |33 



E y (r) 



exp (ik$nsj, 'r 



r± exp(iip) exp (ikons^r 



&0I 



1 — r\ exp(i(p) 



1 - (4vr/3)C" 



,(±) 



sin Q T , 0, ± cos Q?) , f — 2k hn cos 6^, 



(35) 
(36) 
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Figure 1: Coordinate dependences of the amplitude (a, c) and the phase (b, d) of the local 
electric field strength inside the film: - - numerical solution of the system of equations (33), 
■ — ■ — calculation according to the Airy formulae for the case when the film occupies the 
region of space —a(N — 1) < z < 0, — — — calculation according to the Airy formulae for 
the case when the film occupies the region of space —a(N — 1/2) < z < a/2. The calculations 
are performed for the cubic lattice. The external wave is polarized along one of the translation 
vector of the lattice (s-polarization). The parameters are: ko = 0.01 nm" 1 , a = 0.5 nm, N = 20, 
6/ = 30°, C = 0.1 (a, b); C = 0.2 (c, d). 
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where Gt is the angle of refraction, t± and r± are Fresnel coefficients of transmittance and 
reflectance for s-polarization, respectively, h = a(N — 1) is a film thickness, n is a refractive 
index defined by the Lorentz-Lorenz relation (|18|). 

In the latter case h = aN, and in eq.(p5|) one has to make the replacement r — > r — 
na/2, Eqj — > Eqj exp(— ik^a cos 0//2). It seems that the first form of the Airy formulae better 
corresponds to the initial statement of the problem. However, the results of numerical solution 
of the system ( p3|) are in a better agreement with the second form. The reason is that in the 
former case the macroscopic density p which enters the Airy formulae appears to be higher 
than its true value (a\&i x a 2 |)~ 1 , but in the latter case p in the Airy formulae is equal to 
(a | ax x a 2 |) _1 as it should be. 

As it follows from the numerical calculations, the consideration of the discrete structure 
influences the most significantly the behavior of the field amplitude near the film surface where 
rapid oscillations of the local field take place. The amplitude of these oscillations and the 
penetration depth into the film volume increase with the parameter C. One can also see on 
Fig. that the consideration of the discrete structure does not influence much the phase of the 
field. This means that the wave vector of the field remains almost the same at any point inside 
the film. Therefore, the real part of the refractive index of the field is constant. Far from the 
boundaries the distribution of the field is regular and approximately coincides the results given 
by the Airy formulae. 



5 Field on the probe 



N oo 



E(r) = A/(r) - 2mCa £ £ [k pq x (k pq x E° 



exp 



ik pq [r-r J 



When the system of equations (|33| ) is solved and all the quantities E°, j = 1,N are determined, 
one can calculate the field at the probe position: 

„oY 

^exp(zk l r J °). (37) 

j=lp,q= — oo '"PI 

Under the condition flT5| ) the field E(r) equals to the field E p acting on the probe. Let's analyze 
in details the expression 



K 



5.1 Field in the wave zone 

Let's consider first the case when the observation point is in the wave zone, i.e., \z — Zj\ — ► oo. 
In this case the summands with p, q ^ do not contribute to the resultant field due to the 
following reason. All the quantities <j| at p, q ^ are of the order of 2n/a and they are much 
greater than k Q if k^a 1. Therefore, K pq ( p7|) are purely imaginary and all the corresponding 
summands in ( p5p decay rapidly with the increase of the distance from the surface. Only the 
term with koo survives. If z > Zj, koo represents itself the wave vector of the reflected wave kfj, 
otherwise ki 



oo 



E(r) 



E f T ar 



Therefore, in the wave zone we have plane propagating waves: 

| Aj(r) + E£ ar exp (fa R r) , z > 0; 
\ E^ ar exp (zk r) 



-i2nC 



N 



ko cos 0/ 



z < -a(N - 1); 
[k R x (k fl x E°)] exp (-ikirj) 



(38) 



E, 



01 



i2nC- 



k cos 0/ 



N 

E[k x (k xE°)]exp(-2k V;) . 



12 




Figure 2: The difference between the reflected wave field in the near zone, calculated under 
(Ef^ 1 ) and without (E r) the consideration of the discrete structure of the medium, as a func- 
tion of the number of monolayers N which constitute the film. The calculations are performed 
for the cubic lattice. The parameters are: k = 0.01 nm _1 , a = 0.5 nm, Qj = 0°. 



This fields are transverse, because kfjE^ = k E^ ar = 0. The amplitudes E^ ar , E^ ar do not 
depend on the coordinates. 

The discrepancy between the values of the field of the reflected wave calculated under 
consideration of the discrete structure and according to the Airy formulae can be of the order 
of 10%. This discrepancy becomes larger with the decrease of the film thickness and with 
the increase of the parameter C (Fig. |2p. The plots on Fig. lead to the conclusion that at 
high enough film thicknesses (N > 20) the Airy formulae give a satisfactory description of the 
reflected waves in the wave zone. In this situation our method of calculation becomes equivalent 
to that developed in Refs. p5], ^6], ^jj . As for the dependence on the incident angle, the maximal 
discrepancy between the amplitudes of the reflected wave calculated under consideration and 
without consideration of the discrete structure is observed at normal incidence. With the 
increase of the incident angle this discrepancy gradually vanishes (Fig. |3|). In the calculations 
presented on Figs. we used the Airy formulae for the case when the film occupies the region 
of the space —a(N - 1/2) < z < a/2 [pfl 

rj. [1 - exp(i(p)] y 

E m = 2 T- — T ex P (-««oa cos 6/) Eqj. (39) 

1 — r j_ exp{i(p) 



Here we use the same notations as in formula (p5|). Besides, as for the field inside the film, the 
results of our calculations are in a better agreement with the Airy formulae compared to the 
case when the film occupies the region of space —a(N — 1) < z < 0. We would like to note 
also that for p-polarized wave the field of the reflected wave near the Brewster angle does not 



vanish 41 



5.2 Field in the near zone 

If the observation point is near the surface, one can observe a number of characteristic features in 
the field behavior which we are going to discuss now. For the sake of simplicity we shall restrict 
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Figure 3: The dependence of the electric field strength amplitude of the reflected wave E{^ r 
on the incident angle of the external wave Qj in the wave zone for s- (1) and p- (2) polariza- 
tions: - - under consideration of the discrete structure, - - according to the Airy formulae. 
The calculations are performed for the cubic lattice. The parameters are: ko = 0.01 nm^ 1 , 
a = 0.5 nm, N = 3, C = 0.2. 



ourselves, if the opposite is not stated explicitly, by the case of quadratic lattice (a x _L a 2 , 
a x = a 2 = a; g x = g 2 = g = 2-n/a). 



5.2.1 Contributions of propagating and evanescent harmonics 



At a small distance from the surface the summands with p, q ^ in fl37 ) can make a noticeable 



contribution. Let's estimate this contribution considering the following example. As it was 



mentioned above in the discussion of formula (25), it is convenient to single out in eq.(|37|) 
evanescent waves E K = J2 p 2+ q 2 =K A pq exp (^kjjr) with the decay coefficients dependent on k. 
The ratio of the amplitudes of the evanescent harmonics E K for k = (plane reflected wave), 
1, 2, 4 is shown on Fig. |4]. Note that the relative contribution of different E K strongly depends 
on the distance between the probe and the film surface. The behavior shown on Fig. [| is 
typical for the distances greater or equal to one lattice constant. Usually if one takes into 
account the modes with p,q — 0, ±1,±2 the tolerance of the calculations is less than one 
per cent. As it follows from the dependences on Fig. f|, at k^a = tt (a typical value for 
photonic crystals and optical lattices) the monolayers do not fill the evanescent waves radiated 
by the other monolayers. At z = 2a and k^a = 0.005 the ratio of decaying and non-decaying 
field components is about 1% and rapidly decreases with the increase of k a. If a < 1/fco, 
then g ^> ko, and as it follows from eqs.(p5|)-(p7|) the form of the evanescent waves is almost 
independent of the incidence angle of the external wave. Coordinate dependences of the total 
amplitude of the evanescent harmonics is shown on Fig. |5|. 

The dependence of the distance L e , at which the amplitudes of the evanescent and propagat- 
ing components are equal, on the parameter k a is depicted on Fig. |^. Due to the rapid decay 
of the evanescent harmonics the quantity L e can be treated as a distance from the surface up 
to which the contribution of the evanescent component to the resultant field is essential. From 
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Figure 4: The ratio of the evanescent field amplitudes near the surface E K to the amplitude of 
the reflected field in the wave zone E R ar . The calculations are performed for the cubic lattice. 
The parameters are: k = O.Olnm -1 , 6/ = 0, the distance between the observation point and 
the surface is equal to a. 1 - lEjE^l 2 - \E 2 /E f R ar \, 3 - \E 4 /E f R ar \. 
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Figure 5: Coordinate dependences of the ratio of the total amplitude of the evanescent waves 
E ev to the amplitude of the reflected wave field in the wave zone E^ R ar . The calculations are 
performed for the cubic lattice. The parameters are: k = 0.01 nm -1 , a = 0.5 nm, 6/ = 0°. 
The observation point is above an atom of the film. 
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Figure 6: The dependence of the distance L e , at which the amplitudes of the evanescent 
and propagating parts of the reflected wave field are equal to one another, on the parameter 
k a: - - exact numerical calculations, — — — an estimate according to formula (42). The 
calculations are performed for the cubic lattice, 9/ = 0°. 



the dependence shown on Fig. |6] it follows that the contribution of the evanescent component is 
essential at the distances less than two lattice constants. The field radiated by the monolayers, 
which are located at larger distances from the observation point, can be treated as propagating 
plane waves. 

The distance L e can be also analytically estimated on the basis of the following consider- 
ation. As it follows from the dependences depicted on Fig. [| the leading contribution to the 
decaying wave is made by the evanescent wave with k = 1. Therefore, the distance L e can be 
approximately determined from the condition 



\Ei | 



E f R ar 



(40) 



In the case fco« < 1 we have k pq 
([40|) anew as 



from which we get 



->pq 



pql 



Then for a monolayer we can write the condition 



2tt 
2- — exp 

k a 



-2tt— 

a 



h = L^ fin 



In . , , 

2ix \koaJ 



(41) 



(42) 



As it follows from the dependences presented on Fig. |, the estimate ( [4"2"D is in a good agreement 
with the results of exact numerical calculations at small a. 

We can approximately estimate the distance L e for a semi-infinite medium in a similar 
manner assuming that the field of the reflected wave in the wave zone and the field of the 
transmitted wave in the medium are defined by Fresnel reflection and transmission coefficients 
rp and tp, respectively. Then for s-polarized external wave after minor algebra we get 



L e 1 , 

— = — In 

a 2ir 



8vr 2 Ct1 



(l - f c 



(43) 
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In the case of the normal incidence r s F = — (n — l)/(n + 1) and t F = 2/(n + 1). Taking into 
account the Lorentz-Lorenz formula fll8|), we obtain 

— = -*-ln[47r(n + l)]. (44) 

a Z7T 



From eq.(|44j) it folows that if the refractive index n is in the range 1.1 -j- 4.0, then L e /a is in 
the range 0.5 4- 0.7. In fact L e /a has to be higher, because as it follows from the numerical 
calculations in section 4 the local field inside the film in the first atomic layer is higher than 
that in the "bulk" . 

In the case of p-polarized external wave incident on a semi-infinite medium the estimation 
of L e will be given by an equation which has similar structure as ([43"D with t s F and r s F replaced 
by t p F and r F , respectively. In this case near the Brewster angle the reflection can be rather 
small but still noticeable. In such a regime we get the values of L e /a which are somewhat bigger 
than 1. Therefore, our estimations show that the propagating and evanescent harmonics can 
make comparable contributions to the optical response of the medium at the distances from 
the surface of the order of one lattice constant. 



5.2.2 Non-transversality of the total fields of reflected and transmitted waves 

In distinction to the wave zone the polarization vector of the field in the near zone has no a defi- 
nite direction. Its magnitude and the direction depend on all three spatial coordinates. Indeed, 
the polarization vector of a certain evanescent harmonic in the r.h.s of eq.(|37D is perpendicular 
to the wave vector k OT , which in its turn has different orientations depending on p and q. In 
addition, each harmonic has its own decay coefficient K pq and the behavior of the harmonic 
along the surface is defined by its own vector gjj . Carrying out the summation over all the 
harmonics we get the properties of the polarization vector mentioned above. In particular, in 
the case of the normal incidence of the external wave with the polarization vector along the y 
axis the field near the film surface has non- vanishing x- and 2:- components (Fig. [7Ja, c). 

For the purposes of the further analysis let's estimate the ratio of the evanescent com- 
ponent of the reflected wave field to the incident wave field. From the Airy formulae for 
the s-polarization it follows that in the case of ultrathin film the local field inside the film 
is Eq ~ Eqj/[1 — (47r/3)C]. The evanescent field at the distance of one lattice constant 
from the surface is mainly determined by the harmonic E\. For the square lattice E\ = 
-2(2tt) 2 C exp(-27r)£ « 0. 147(7/ [1 - (4n/3)C]E 0I . Even for C = 0.2 the ratio E QI « 0.2. 
Thus, the amplitude of the evanescent wave at the distance from the surface of the order of 
one lattice constant is usually much less than the incident field amplitude. Therefore, one can 
conclude that in the field of the probe the main contribution is made by the component of the 
field directed along the polarization vector of the external wave. In particular, comparing the 
intensities of x and ^-components of the field (Fig. 0a, c) with y-component (Fig. 0b) one can 
notice that the contribution of x- and z-components is small and the field behavior near the 
surface is mainly determined by ^-component. 

One can come to one more conclusion regarding the direction of the field polarization vector 
near the film surface. Under the normal incidence of the external wave in the case when there 
are no parallel shifts of the atomic planes relative to one another, i.e., ay = a2j = 0, j = 1, N, 
at the observation points above the atoms the ^-component of the reflected wave field vanishes 
(Fig. 0). Indeed, taking into account the properties of the tensor / and summing in fl3"3"|) the 
harmonics with equal K pq {n pq = K- p - q ), we get that (E°) X) j, are independent of (E°)^. The z 
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Figure 7: The ratio of the intensities of x- (a), y- (b) and z- (c) components of the electro- 
magnetic field on the probe (I x ~ \E X \ 2 , I y ~ \E y \ 2 , I z ~ \E Z \ 2 ) to the intensity of the incident 
wave Io on the coordinates x and y. The parameters are: the number of monolayers in the film 
iV = 20, C = 0.2, k = 0.01 nm" 1 , a = 0.5 nm, d = a, 0/ = 0°, $/ = 0°. The intensity minima 
correspond to the atomic positions. In spite of the fact that the external wave is polarized 
along the y axis, optical responds of the medium contains all three components. 
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components of the field inside the film satisfy to the system of homogeneous equations which 
has only the trivial solution (E°) 2 = 0, j — 1,N. Then it follows immediately from (jS?D that 
at the observation points above the atoms (E p ) z = 0. 

5.2.3 Dependence of field intensity on the longitudinal coordinates x and y 

The electromagnetic field intensity is defined by its strength squared. In a general case it can be 
calculated numerically. But in the case of the normal incidence (G/ = 0) one can obtain rather 
simple analytical expressions. Let's consider this case. From the plots on Fig. |6] it follows 
that the component of the field, which is a periodic function of longitudinal coordinates, is 
induced mainly by the surface atoms. The bulk atoms produce only a constant background. 
Thus, in order to analyze the behavior of the evanescent waves it is enough to consider only 
one surface monolayer (zj = in eq . (|37|) ) . For the qualitative analysis of the field behavior 
near the film surface it is enough to keep in eq.fl3"?D only the harmonics with (p,q) = (1,0), 
(— 1, 0), (0, 1), (0, —1) (we assume that the quantities g\ and <?2 differ from one another no more 
that two times). The amplitudes of all other harmonics are at least one order of magnitude 
less [ 42| . Omitting the unimportant constant factor after an elementary algebra we obtain an 



approximate expression for the field intensity at the probe position: 

/ « L/ + Re(4 ar ) (45) 



_ ^ C a £ {Re( gi E?) cos [ B (r - r°)] - g t Re(E° lz ) sin [ a (r - r°)] } 



=1,2 E 0I 



expl-gjZj 
9i 



1 2 



which is valid for the films of arbitrary thickness. Here E® is an amplitude of the unperturbed 
part of the field in the first (surface) monolayer. In fact, this expression is the second power 
of the real part of the field component at the probe position, directed along the polarization 
vector of the external field. The imaginary part of this component and all other components 
are negligible. 

From the expression (HS[) one can see that in general the intensity minima and maxima 
are shifted relative to the atomic positions. This is due to the term which contains E\ z . Due 
to the fact that even in the case of the normal incidence the dependence of the direction 
of the vector on the direction of Eo/ is rather complicated, the location of the intensity 
minima and maxima depends on E / in a complicated manner as well. Let's consider the 
case when there are no parallel shifts of the atomic planes relative to one another. As it was 
shown above, in this case E\ z = and the expression ([45"D becomes simpler. In this case the 
intensity minima are located exactly at the positions of the atoms. Indeed, at these points 
cos[gi(r — 1*5)] = cos[g2(r — r±)] = 1. Thus, formula fl45|) explains the contrast reversal, which 



has been obtained in Refs. [^, [E| by means of numerical calculations. 

In analogous manner one could assume that the intensity maxima are located in the centers 
of the elementary cells, where the cosines are equal to —1. However, this is not exactly the 
case as one can see on Fig. [8]. With the change of the mutual orientation of the polarization 
vector of the external wave and the lattice elementary translations vectors the location of the 
intensity maxima changes as well. The reason is the following. Let the polarization vector of 
the external wave is perpendicular to one of the vectors of the reciprocal lattice, for instance, 
(E /gi) = 0. The field behavior along the direction &i in this case is determined by higher 
order harmonics with the smallest decay coefficient and with non-vanishing component along 
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Figure 8: The ratio of the electromagnetic field intensity on the probe / to the incident wave 
intensity 1$. The parameters are: N = 20, C = 0.2, /c = 0.01 nm" 1 , a = 0.5 nm, d = a, 
0/ = 0°; $/ = 0° (a), 30° (b), 45° (c). The electromagnetic field intensity is a periodic function 
of the coordinates x and y, and the period is equal to the lattice constant a. The intensity 
minima correspond to the atomic positions. 
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the polarization vector of the external wave. Obviously, these are harmonics with \p\ 
In eq.(H^) in the square brackets one has to add one more term 
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where g™ n is one of the vectors g pq with \p\ = \q\ = 1, which has the minimal length. At the 
observation points above the atoms (r" = rj[ ) this term is minimal, i.e., the intensity minima 
are still above the atoms, even provided that (E /gi) = or (E^gi) = 0. The locations of the 
intensity maxima can be found from the condition (g™ n r) = n ± 2irk, where k is an integer. 
This condition fulfills not in the centers of the elementary cells, like in the case when the 
directions of the polarization vector of the external wave Eo/ and the local field inside the film 
E^ does not coincide with any of the translation vectors of the lattice (Fig. ||b,c), but at the 
points r" = rj[ +aj/2 or r" = rj[ +a2/2 (Fig. |8]a). Thus, in distinction to the intensity minima 
locations, which are always above the atoms, the locations of the intensity maxima depend on 
the mutual orientation of the field polarization and the translation vectors of the lattice. 

Formula ([45]) allows one to investigate the image contrast, which is defined here as a maximal 



difference of the light intensity along a certain direction. In Refs. j|, [36], [44| it was point 
out that the field component parallel to the scan direction displays a better contrast than 
that perpendicular to this direction. Let's consider this problem in more details. The image 
contrast in the direction aj is determined by the mutual orientation of the vectors Eo/, E°, and 
gj. The maximal contrast is reached at some angle between E / and gj when the coefficient 
(gjE /)Re(g i E5 ) ) takes the largest value. This angle depends in general on the symmetry of the 
atomic distribution at the surface and on the atomic polarizability. We are going to come back 
to this issue in our subsequent publications. 

The minimal contrast along the direction ai for the lines which pass through the atoms 
is achieved when the polarization vector of the external wave or the field vector in the film 
is perpendicular to the vector of the reciprocal lattice gi (parallel to the vector of the real 
lattice a 2 ) (f45"l). Performing numerical calculations, one can show that in general the conditions 
(gjEo/) = and (gjE°) = are satisfied at different directions of the polarization vector of 
the external wave. Thus, in general there exist two angles at which contrast minima can be 
observed. The angle between the two directions of the polarization vector of the external field 
allows one to get the information about the local field inside the film and, therefore, about the 
atomic polarizability. 

Another situation occurs for the contrast along the lines which pass through the centers 
of the elementary cells. In this case the contrast minima are determined by the condition of 
mutual compensation of the contributions from the harmonics with decay coefficients gi and 
5™. The compensation occurs at some angle <p c between the polarization vector of the field 
inside the film and the vector g 1; which depends on the symmetry of the atomic distribution 
at the surface. The angle 4> c strongly depends on the distance from the surface. 

Let's consider the simplest case of the cubic lattice. Under the normal incidence of the 
external wave the direction of the polarization vector of the field inside the film coincides with 
that of the external field. In this case there are two directions of the polarization vector of 
the external wave, for which one can observe the contrast minima along the lines which pass 
through the atoms. This corresponds to the situations when the polarization vector coincides 
with ai or a 2 . For the cubic lattice one can get an explicit expression for the angle C at which 
one observes the minimal contrast along the lines which pass through the centers of the cells. 
The dependence of the angle <p c on the distance between the observation point and the surface 
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Figure 9: The angle between the minimal contrast direction and the polarization vector of 
the external wave <p c as a function of the distance to the film surface d. The calculations are 
performed for the cubic lattice, 0/ = 0°. 



is given by the approximate formula 

sin C = 2 1/4 exp [-tt(V2 - l)d/a] , (46) 

which is valid only for d > a. The dependence 4> c {d) for the cubic lattice is shown on Fig. |^. As 
one can see on the figure, the angle <p c can take the values of dozens of degrees at the distance 
from the surface of the order of one lattice constant. The magnitude of the angle remains 
essential at the distances of several lattice constants. 

Thus, changing the direction of the polarization vector of the external wave one can ma- 
nipulate the image contrast. One can determine the atomic positions and the symmetry of the 
crystal lattice detecting the intensity minima. Measuring the angle between the directions of 
the polarization vector of the external wave, at which the minimal image contrast is reached, 
one can get the information about the local fields inside the film. Having on hand this in- 
formation and measuring the angle C between the translation vectors of the lattice and the 
polarization vector of the external wave one can determine the distance between the probe and 
the surface. 

The analysis based on eq.(p5[) allows one to say that formula ( f45|) remains valid in the case 



of arbitrary incident angle of s-polarized wave as well. One has to note that in the near-field 



optical measurements it is more useful to work with p-polarized waves | 36fl . In the case of 
arbitrary incidence angle of p-polarized wave the distribution of the field in the near zone is 
essentially different. This has been shown in Ref. || on the basis of numerical calculations. We 
intend to give a detailed description of the field behavior in the near zone at arbitrary incident 
angles in our subsequent publications. 



6 Conlcusions 

We have solved a boundary problem of the linear classical optics devoted to the investigation of 
the electromagnetic field behavior near the surface of a dielectric medium taking into account 
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its discrete structure. The main attention has been paid to the investigation of the near-zone 
optical properties of dielectrics. It has been shown that at the distances from the surface 
less than two lattice constants the behavior of the reflected and transmitted waves is entirely 
different from what we have in the wave zone. The intensity distribution in the near zone allows 
one to determine the atomic positions at the surface. 

We have shown that at the distances from the surface larger than the interatomic distance in 
the film the probe which measures the field does not significantly influence the field distribution 
in the film. In the present paper we have suggested a method of the optical control of the 
distance between the probe and the surface. The method is based on the analysis of the near- 
field distribution at different directions of the polarization vector of the external radiation. 

At first glance the system treated in our paper (ultrathin dielectric film without any sub- 
strate) seems to be unrealistic. However, it can be used for the investigation of the near-zone 
optical properties of the films of arbitrary thicknesses. As we have shown the near- zone optical 
response of a material is defined mainly by the evanescent harmonics, which play an important 
role only within few interatomic distances near the surface, because they decay very rapidly. 
This means that the "bulk" monolayers influence significantly only the wave-zone optical prop- 
erties and does not change much the ratio between propagating and evanescent harmonics. 

Our results can be useful not only for the development of the ultrahigh resolution near-field 
microscopy, but also for the investigation of the optical properties of photonic crystals in the 
long wavelength approximation. 

The results obtained in the present paper are valid not only in the case of an ideal crystal 
lattice, but also when the medium has only short-range order as in the case of real surfaces. 
Indeed, the methods based on the Fourier transform and the Lorentz method give analogous 
results. At the same time only the short-range order is important for the Lorentz method and 
the magnitude of the dipole field is mainly determined by the atoms located at the distance 
from the observation point of the order of 2 4- 3 lattice constants. The atomic distribution 
outside of this region does not practically influence the field. 
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